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1 Introduction 
1.1 Overview 

In this paper we study surfaces in that arise as limit shapes in random sur- 
face models related to planar dimers. These limit shapes are surface tension 
minimizers, that is, they minimize a functional of the form J a{Vh) dx dy 
among all Lipschitz functions h taking given values on the boundary of the 
domain. The surface tension a has singularities and is not strictly convex, 
which leads to formation of facets and edges in the limit shapes. 

We find a change of variables that reduces the Euler-Lagrange equation for 
the variational problem to the complex inviscid Burgers equation (complex 
Hopf equation). The equation can thus be solved in terms of an arbitrary 
holomorphic function, which is somewhat similar in spirit to Weierstrass 
parametrization of minimal surfaces. We further show that for a natural 
dense set of boundary conditions, the holomorphic function in question is. 
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in fact, algebraic. The tools of algebraic geometry can thus be brought in 
to study the minimizers and, especially, the formation of their singularities. 
This is illustrated by several explicitly computed examples. 

1.2 Random surface models 

The general class of random surface models to which our results apply are 
the height functions of dimer models on periodic weighted planar bipartite 
graphs, see for an introduction. For our purposes, the following simplest 
model in this class is already very interesting. 

Consider continuous surfaces in glued out of unit squares in coordinate 
directions, that is, out of parallel translates of the faces of the unit cube. 
We assume that the surfaces are monotone, that is, project 1-to-l in the 
(1,1,1) direction, and span a given polygonal contour C, see Figured] for an 
illustration. We call such surfaces stepped surfaces. Uniform measure on all 
such stepped surfaces is an obvious 2-dimensional generalization of the simple 
random walk on Z. It is well-known and obvious from Figure [T] that stepped 
surfaces spanning a given contour C are in a bijection with tilings of the 
region VL enclosed by the (1, 1, l)-projection of C with three kinds of rhombi, 
known as lozenges. They are also in bijection with dimer coverings — that is, 
perfect matchings — of the hexagonal graph. 

Let Cn be a sequence of boundary contours each of which can be spanned 
by at least one stepped surface; we will call such contours tilahle. Suppose 
that Cn converges to a given curve C C M?. In this case C can be spanned 
by a Lipschitz surface whose normal (which exists almost everywhere) points 
into the positive octant (any limit of stepped surfaces of the C„ is such a 
surface); we will call such contours C feasible. It is well-known, see e.g. [H] 
that feasibility is equivalent to being a limit of tilable contours. 

It is natural to consider the weak limit of the corresponding uniform 
measures on stepped surfaces of C„, scaled by n~^. A form of the law of 
large numbers, proved in j3j , says that this limit depends only on the contour 
C and is, in fact, a 5-measure on a single Lipschitz surface spanning C. 
This surface, known as the limit shape, can, for example, be interpreted as 
the macroscopic shape of an interface (e.g. crystal surface) arising from its 
microscopic description and boundary conditions. 

A question of obvious interest is to describe this limit shape quantita- 
tively and qualitatively. In particular, an important phenomenon, observed 
in nature and reproduced in this simple model, is formation of facets, that 
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Figure 1: A limit shape simulation. The "frozen boundary" (yellow), which 
is asymptotically the boundary of the facets, is a cardioid. 
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is, flat pieces, in the limit shape. 



1.3 Variational principle 

The following variational characterization of the limit shape was proved in 
[3]. The limit shape may be parameterized as the graph of a function 

X3 = h{x3 - Xi, X3 - X2) 

where /i is a Lipschitz function with gradient Vh (which is defined almost 
everywhere by Rademacher's theorem) in the triangle 

A = conv{(0,0),(0,l),(l,0)}. (1) 

Let fl be the plane region enclosed by the projection of the boundary curve 
C along the (1, 1, l)-direction. We will use 

{X,y) = (X3 - Xi,X3 - X2) 

as coordinates on Q. The limit shape height function h is the unique mini- 
mizer of the functional 

ENT(/i) = / a{Vh{x,y))dxdy ^ mm (2) 
Jn 

among all Lipschitz functions with the given boundary conditions (that is, 
for which the graph of h, restricted to dfl, is C) and gradient in A. Here a 
is the surface tension. It has an explicit form which, in the language of |13] , 
is the Legendre dual of the Ronkin function of the simplest plane curve 

z + w = 1 . 

We recall that for a plane curve P{z,w) = 0, its Ronkin function [TJ] is 
defined by 

[2tXiY J J\z\=e^ I ' Z W 

The gradient always takes values in the Newton polygon A(P) of the 
polynomial P, so A(P) is naturally the domain of the Legendre transform 
a = "JV^ . For the straight line as above, the Newton polygon is evidently the 
triangle ([T]). See Figured 
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Figure 2: Graph of —a for the hne w) = z + w — 1. 
1.4 General planar graphs 

More generally, to an arbitrary planar periodic bipartite edge-weighted graph 
(with periodic edge weights) one associates a real plane curve P{z,w) = 0, 
its spectral curve. Limit shapes for the height function in the corresponding 
dimer problem satisfy the same variational principle as above with a = 3^^, 
where 3^ is the Ronkin function of P{z,w), see [13]. The spectral curve P 
of the dimer model is always a so-called Harnack curve. This property, also 
known as maximality, is crucial for the constructions in this paper. Har- 
nack curves form a very special class of real plane curves, see [1^. Many 
probabilistic implications of maximality are discussed in |13] . 

Maximality, in particular, implies that the surface tension cr is analytic 
and strictly convex in the interior A(P)° of A(P), with the exception of 
lattice points, that is, points of Z^. At a lattice point, generically, a has 
a conical singularity, often referred to as a cusp in physics literature. For 
example, at the point (0, 0) in Figure [21 a has the form 

t s 

a(s, t) = tin h s In h o(s, t) , 

^ ^ t+s t+s ^ ^ 

in which the reader will recognize the Shannon entropy. More generally, 
at a lattice point, the graph of a is approximated by a cone over a curve 
which is the polar dual of the curve bounding the corresponding facet of the 
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Ronkin function. Such singularities lead to facet formation in the minimizer. 
They disappear only if the spectral curve has a nodal singularity, which is a 
codimension two condition on the weights in the dimer model. 

At the boundary of the Newton polygon, a is piecewise-linear and, in 
particular, not strictly convex. This loss of strict convexity leads to edges in 
the limit shape. Both facets and edges are clearly visible in the simulation 
in Figure m 

These features of the surface tension make the analytic study of the vari- 
ational problem ([2]) challenging and interesting. 



1.5 Complex Burgers equation 

We say that a point (a;o, yo) G is in the liquid region if the minimizer h{x, y) 
is in a neighborhood of (xq, yo) and Vh{xQ, yo) is either in A(P)°\Z^ or at 
a point (in Z^) corresponding to a nodal singularity of P. This terminology 
is motivated by the frozen/liquid/gaseous classification of phases of a dimer 
model [T3] . 

Since a is analytic and strictly convex on the liquid region, Morrey's 
theorem [T7] implies that h is analytic there and satisfies the Euler-Lagrange 
equation div Vcr o V/i = 0. A more general equation 

divVdo V/i = c, (3) 

where c is a real constant, arises in the volume- constrained minimization 
problem 

= / <y(yh{x,y))dxdy + c I h{x,y) dx dy ^ mm , (4) 
Jn Jn 

in which c is the Lagrange multiplier for the constraint 

h{x,y) dx dy = const 

that forces the limit shape to enclose a given volume (or more precisely, to 
contain a given volume under its graph). 

Our first result is the following reduction of the Euler-Lagrange equation 
(I3D to a first-order quasilinear equation for a complex valued function z{x, y). 
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Theorem 1. In the liquid region, we have 



V/i = — (argw, — argz) (5) 

TT 



where the functions z and w solve the differential equation 



^ + ^ = c (6) 

z w 

and the algebraic equation P{z,w) = . 

Note that in ([5]) the argument of a complex number is a muhivalued 
function and the statement is that there exists a branch of the argument for 
which equahty ([5]) is satisfied. 

For example, in the lozenge tiling case 

P{z,w) = z + w -1 (7) 
and without the volume constraint, equation ([6]) becomes 



Zy 1 — z 

which becomes the standard complex inviscid Burgers equation 

Ua: = UUy (8) 

upon the substitution u = zjil — z). In general, equation ([6]), being a first- 
order quasilinear equation for z{x,y), has properties parallel to those of ([8]) 
and, in particular, it can be solved by complex characteristics as follows: 

Corollary 1. For c ^ there exists an analytic function Q of two variables 
such that 

Q{e-'^z,e-'yw) = 

in the liquid region. When c = there is an analytic function Qo of two 
variables such that 

Qo{z, w) = xzP^ + ywP^. (9) 

In other words (for c 7^ 0), as functions of x and y, z and w are found by 
solving the system 

= 0. 

Q z, e-^y w) = 0. ^ ' 
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A general solution to the Euler-Lagrange equation is thus parameterized 
by an analytic curve Q = in the plane. 

Below we define a dense set of boundary conditions for which Q is in fact 
algebraic. In this case ffTOj) is a system of algebraic equations and, hence, h 
is an integral of an algebraic function of e~'^^ and e~^^. 

1.6 Frozen boundary 

The liquid region may extend all the way to the boundary of Q, but it may 
happen that it ends at a facet, that is, a domain where the limit shape height 
function h is a linear function with integral gradient, as in Figure [H We call 
the boundary between the liquid region and the facets the frozen boundary. 

There are two kinds of facets: the ones with gradient V/i in the interior of 
the Newton polygon A(P) and the ones with V/i on the boundary of A(P). 
While the corresponding microscopic properties of random surfaces are, in a 
sense, diametrically opposite (they correspond to gaseous and frozen phases, 
respectively, in the classification of [13]), the only difference relevant for this 
paper is that only gaseous facets can occur in the interior of the domain. We 
will call such floating facets bubbles. 

In terms of ([5]), a point {z, w) of the frozen boundary is a real point of the 
spectral curve P{z, w) = 0. A curve's worth of such real solutions of (fTOj) is 
possible only if Q is real. In this case, at the frozen boundary, the solutions 
{z, w) and {z, w) of the system (ITU]) become identical. That is, the frozen 
boundary is a shock for the equation ([6]). In contrast to the usual real- valued 
Burgers equation, a shock is a special event in the complex case as complex 
characteristics tend to miss each other. 

Given two analytic functions P and Q, the locus where the system f llUI) 
has a double root is itself an analytic curve in the (a;, y)-plane. In the case 
when P and Q are algebraic, this discriminant curve is of the form 

^(e-^^e-'^J') = 0, (11) 

where R is polynomial which can be written down effectively in terms of 
resultants. 

The frozen boundary i? is a subset (in principle, proper) of the discrimi- 
nant curve fill I) . A more invariant way of saying the same thing is that the 
map 

£xp : (xi,X2,a;3) ^ (e^^^^ : e"^^ : e"^^) G MP^ 
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from to the real projective plane takes the frozen boundary to a real 
algebraic curve. As c — the map £xp, suitably rescaled, becomes the 
projection in the (1, 1, 1) direction. Monotonicity insures that our random 
surfaces are mapped 1:1 by £xp. Indeed, if (xi, X2, x^) is on the surface then 
(xi +t,X2 + t, X3 + t) is on the surface only for t = 0. 

At a generic point, the frozen boundary is smooth and flTU]) has exactly 
one real double root. The functions z{x,y) and w{x,y) have a square-root 
singularity at the frozen boundary and, hence, the arguments of z and w grow 
like the square root of the distance as one moves from the frozen boundary 
inside the liquid region. Integrating ([5]), we conclude that the limit shape 
height function has an x^/^ singularity at the frozen boundary. Thus one 
recovers the well-known Pokrovsky-Talapov law [2T] in this situation. 

At special points of the frozen boundary, triple solutions of ffTOl) occur. 
At such points, the frozen boundary has a cusp singularity. One such point 
can be seen in Figure [TJ 

1.7 Algebraic solutions 

Any boundary contour can be approximated by piecewise linear contours and 
even by piecewise linear contours C with segments in coordinate directions. 
For such contours the function Q is algebraic (of degree growing with the 
number of segments). In this paper we will prove this in the simplest case 
when P{z, w) = z + w — 1 and the boundary C is connected. The precise 
result that we will prove is the following 

Theorem 2. Suppose P{z,w) = z + w — 1 and the boundary contour C 
is feasible, connected, and polygonal with 3d sides in coordinate directions 
(cyclically repeated). Then Q is an algebraic curve of degree d and genus 
zero. 

The corresponding result for cases in which some of the edge lengths are 
allowed to be zero can be obtained as limits, letting the edge lengths shrink. 
In these cases Q may have smaller degree. 

1.8 Reconstruction of limit shape from the Burgers 
equation 

The frozen boundary i? is a genus zero curve inscribed in the polygon £xp(C). 
The reconstruction of the limit shape from R has the following elementary 
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geometric interpretation. 

Let X = {xi,X2,X3) be a point on the limit shape in the hquid region. 
Then £xp(x) hes inside the frozen boundary R. In the course of proving 
Theorem [2l we will show that there is a unique, up to conjugation, complex 
tangent to R through the point £xp(x). This is illustrated in the Figure [31 
There are 3 tangent lines to a cardioid through a general point in the plane 
(this number of tangents is the degree of the dual curve and is classically 
known as the class of a curve). Through a point outside the cardioid, all 
3 of these lines are real. Through a point inside the cardioid, only one real 
tangent exists, the other two tangents are complex. 




Figure 3: Tangent lines to a cardioid 
Let the complex tangent through £xp(x) be given by the equation 

aiZi + 02-22 + 0-3-23 = 

in homogeneous coordinates {zi : Z2 '■ z^) on the projective plane. Since the 
point £xp(x) lies on this line, this gives us a triangle in the complex plane 
illustrated in Figure HI 




Figure 4: A complex tangent defines a triangle 
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The numbers (01,02,03) are unique up to common complex factor and 
complex conjugation. If follows that the similarity class of the triangle in 
Figure mis well defined, in particular, its angles 01,02,0:3 are well defined. 
Formula ([5]) simply says that 

(01,02,03) e ]R>o 

is the normal to the limit shape at the point x. 

As the point x approaches the frozen boundary, the complex tangent 
approaches the corresponding real tangent, the triangle degenerates, and the 
normal starts pointing in one of the coordinate directions. 

Higher genus frozen boundaries occur for multiply connected liquid re- 
gions. The holes may be created manually by considering surfaces spanning 
a disconnected boundary contour or they may appear as bubbles for spectral 
curves P of genus > 1. Examples of such situations will be considered in 
Sections 15.11 and 15. 2[ In the multiply-connected case, the Burgers equation 
needs to be supplemented by period conditions, see Sections 12. 5[ 

1.9 Burgers equation in random matrix theory 

In the case when Q is an infinite vertical strip, the random surface model 
studied here can be conveniently restated in terms of nonintersecting lat- 
tice paths connecting given points on dQ. The continuous limit of such 
nonintersecting path models is well-known to be described by the Hermitian 
2-matrix model (in other words, by a random walk on the space of Hermitian 
matrices with boundary conditions in given conjugacy classes). In physics 
literature, complex Burgers equation first appeared in this context in the 
work of Matytsin [H]. Its rigorous mathematical study appears in the work 
of A. Guionnet [5]. Generalizations of these equations to interacting particle 
systems are discussed in [1]. 
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2 Complex Burgers equation 
2.1 Proof of Theorem [1] 

By tlie basic properties of the Ronkin function, see [T5], the function 

V(T = (va^)~^ : ^ R2 

takes values in the amoeba of the spectral curve. By definition, this means 
that for any {x,y) in the hquid region we can find a point {z,w) on the 
spectral curve satisfying 

[Va o V/i] (x, y) = (log 1^1, log \w\) . (12) 

By our hypothesis, near (x, y) the LHS of (fT2|) lies in the interior of the 
amoeba. One of the several equivalent definitions of a Harnack curve is that 
the map 

(z, w) (— > (log |z|, log |w|) (13) 

from the curve to its amoeba is 2-to-l with nonzero Jacobian over the amoeba's 
interior. Therefore, the lift (x, y) i-^ {z, w) in (fT2|) is unique up to complex 
conjugation, and real analytic in a neighborhood of {xo,yo)- 

Using (fT2|) . the Euler-Lagrange equation ([3]) stated in terms of {z,w) 
becomes 

\ z w / 

This first order equation for the gradient V/i needs to be supplemented by 
the usual consistency relation for the mixed partials: 

Again, it is a special property of the Harnack curve that 

VJi (log \z\, log \w\) = — (t arg w, ± arg^;) (15) 

TT 

where on the right we mean a lift of (=F argt;;, ± arg^;) and where the choice 
of signs depends on the choice of one of the preimages in f|T3|) . see [20] and 
also [T5]. This turns the relation f|T^ into the imaginary part of ([6]) (since 
c G M) and concludes the proof. 
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2.2 Proof of Corollary [T] 

Since z and w satisfy an algebraic equation P{z, w) = 0, their Jacobian 
vanishes 

J{z,w) = det(^'' ^A=0. 

We have 

J{e~^^z^ e~^'^w) — J{z, w) = c^zw — czwy — cz^w = 

by and hence J{e~'^^z, e~'^^w) = identically. The functions e~'^^z and 
e~^^iy are real analytic in the liquid region and the vanishing of the Jacobian 
implies a functional dependence of the form 

g(e-^^2,e-^^w) = 0, 

where Q is an analytic function (note that if ly = W{x,y),Z = Z{x,y) 
satisfy Z^Wy = ZyW^, then dW/dZ = W^/Z^ - Wy/Zy = 0, so that W is an 
analytic function of Z, and vice versa.) 

2.3 Complex structure on the liquid region 

Consider the function 

C = e-'^^z 

which is defined on the liquid region. Differentiating P{z,w) with respect to 
y and using ([6]) implies that 

when P{z, w) = 0. The rational function 

{z,w)^j{z) = ^eCP' 

is known as the logarithmic Gaufi map of the spectral curve P{z,w) = 0. 

By a result of Mikhalkin [15], the preimage 7"^ (MP^) is precisely the 
critical locus of the map from the curve to its amoeba. For a Harnack curve, 
therefore, it coincides with the real locus. It follows that on the liquid region 
7 takes values in upper or lower half-plane, depending on the choice of lift of 
( IT2|) . With one of the choices, the map 

{x,y) C 
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is orientation-preserving and we can use it to pull back the complex structure. 

Note that one could equally well take the function e~'^'^w as defining the 
complex structure on the liquid region. However, by Corollary [1] the complex 
structure thus obtained is the same. Yet another way to say this is that 

maps the liquid region to the analytic curve Q = and the complex structure 
on the liquid region is the unique one making this map holomorphic. 

The orientation-preserving property of Ci^jV) can be strengthened as 
follows 

Proposition 2. For c ^ and {x,y) in the liquid region, the map 

nQ:ix,y)^ie~'-z,e-'yw)eQ (17) 

is an orientation preserving local diffeomorphism onto its image. In the case 
P{z,w) = z + w — 1 the map is a diffeomorphism onto its image. 

Proof. Suppose that ttq maps two distinct points {xi,yi) and {x2,y2) to the 
same point 

of Q, where Zi = z{xi, yi). This means that on the curve P we have a pair of 
distinct points 

P 3 ^ {Z2,W2) e P 

with the same arguments 

axg{zi,wi) = aTg{z2,W2) . 

However, P is a Harnack curve and for a such a curve we have (ITSll . By 
strict convexity of the Ronkin function, its gradient maps both components 
of -P(C) \ P(R) one-to-one to the Newton polygon of P, and, as a result, a 
point of P is locally determined by the arguments. l{P = z + w — 1 (or 
even if P has degree < 2), a point of P in a component of P(C) \ P(M) is 
determined uniquely by its arguments. □ 

There are examples where ttq is a non-trivial covering map. 

In the case c = 0, we can use 2; or to define the complex structure. 
This complex structure is used in |T0] in the study of the fiuctuations of the 
random surfaces. Conjecturally, the Gaussian correction to the limit shape 
is given by the massless free field in the conformal structure defined by (. 
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2.4 Genus of Q 



Suppose that the hquid region is surrounded by the frozen boundary, that 
is, suppose that Qz, — at the boundary of the hquid region. The real 
parts and may or may not be continuous at the boundary of the 
hquid region. 

Proposition 3. If the liquid region is surrounded by the frozen boundary and 
the functions z{x,y) and w{x,y) are continuous up to the frozen boundary 
then the curve Q is algebraic. IfP = z + w — 1 then Q is of the same genus 
as the liquid region. 

Proof. If z,w are continuous then the map f[T7l) . and its complex conjugate 
glue along the boundary of the liquid region to a map from the double of 
the liquid region to the curve Q. This map is orientation preserving and 
unramified by Proposition [2], and so a covering map of Q. It follows that the 
analytic curve Q is compact. lfP = z + w — l, the map is a diffeomorphism 
and therefore Q is of the same genus as the liquid region. □ 

Conversely, if Q is algebraic then the solutions of (ITO!) are continuous. 
However, only for some very special curves Q will the solution of (fTOj) satisfy 
the global requirements of being single- valued in some region and becoming 
real at the region's boundary. 

By taking limits of algebraic solutions, one can construct solutions having, 
for example, accumulation points of cusps. At such a point z is discontinuous. 
Namely, z turns infinitely many times around MP^ as one approaches the 
accumulation point along the frozen boundary. 

It would be interesting to prove the continuity required in Proposition [3] 
directly from the variational problem and polygonality of the boundary con- 
ditions. In this paper, we take a different route and construct the algebraic 
curve Q by a deformation argument. 

Note that Proposition [3] forces Q to have the maximal possible number 
of real ovals (real components), namely genus plus one. Such real curves are 
known as M-curves. 

2.5 Islands and bubbles 

Holes in the liquid region may arise by considering surfaces spanning a dis- 
connected boundary contour or they may appear as gaseous bubbles (where 
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the surface tension has a cusp). We will call such holes islands and bubbles, 
respectively. A schematic representation of an island and a bubble can be 
seen in Figure [5] 



Figure 5: An island, a bubble, and their period conditions. 

In principle, the solutions of a system like (ITOl) may pick up a monodromy 
around a nontrivial loop in the liquid region. However, in our case we require 
them to be single-valued because they are determined by the single-valued 
function V/i. 

In the presence of holes, the Burgers equation ([6]) needs to be supple- 
ment by integral constraints. This constraints can be phrased as imposing 
conditions on the periods of the curve Q. 

2.5.1 Areas of bubbles 

In the case of a bubble inside the liquid region, it is more convenient to use 
the integrated form of ([3]) which says that the flow of the vector field Vcro V/i 
through any closed contour equals c times the area enclosed. 

Let a domain B G Q he such that its boundary dB is smooth and lies in 
the liquid region. Consider a variation of h which is a smooth approximation 
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to the characteristic function of B. This yields the real part of the equation 

/ a = c Area(5) . (18) 

JdB 

where a is the following 1-form 

a = log zdy — log w dx . 

By the vanishing of the imaginary part of flTSl) means that the height 
function h is well-defined. Thus flTK]) is satisfied by any B as above. 

Suppose that {z, w) are continuous up to the frozen boundary. On the 
frozen boundary, the points {e~'^^z, e~'^^w) and {z, w) are real points of Q and 
P, respectively. This allows us to associate, to each oval of Q, a corresponding 
oval of P. We will call an oval compact if it does not intersect the coordinate 
axes. Note that an oval of Q is compact if and only if the corresponding oval 
of P is. Because P is a Harnack curve, it has exactly one noncompact oval 
that intersects each coordinate axis degP times. 

If P is a bubble, a point on its boundary gives a point {z, w) & P belonging 
to a compact oval of P. Integrating f[T5]) around the boundary of B, the 
imaginary part of the left-hand side of f[T5]) is then zero, since on the boundary 
of a bubble z and w are real and of definite sign. However the real part 
may not vanish and represents a nontrivial constraint of the curve Q. This 
constraint can be reformulated as follows. 

Given a compact oval O of an plane curve, consider its logarithmic area 

LogArea(O) = / \ogzd\ogw, 
Jo 

that is, the area enclosed by the image of O under the amoeba map. 

Proposition 4. The logarithmic areas of compact ovals of Q are multiples 
of the logarithmic areas of the corresponding ovals of P. The multiple at a 
given oval is given by the winding number of the map from the corresponding 
bubble B to the oval of Q. 

Proof. Let Op and Oq be the ovals of P and Q corresponding to a bubble 
B. Because these ovals are compact, the logarithms of z, w, ( = e~^^z and 
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uj = e '^^w are well-defined on dB. We have 
LogArea(OQ) = / logCc?logu; = 

JdB 

xdy — c {y dlogz + X dlogw) + / logz dlogw = 

JdB JdB JdB 

(? Area(i?) — c / a + m LogArea(Op) = m LogArea(Op) , 

JdB 

where m is the winding number of dB around Oq. □ 
2.5.2 Heights of islands 

Different period conditions need to be imposed if VL has interior boundary as 
in Figure or in the example considered in Section [F!T1 In this case, z and w 
take values on the big oval of the spectral curve and do change sign on dB. 
In this case, the vanishing of the imaginary part of (fT8l) is equivalent to the 
condition that the polygon inscribed in B lifts to a closed contour in M^, see 
the discussion after the statement of Theorem [3] below. The argument that 
was used to derive the real part of ( fTSl) is no longer valid and, in fact, the real 
part of (ITSll no longer holds. Instead, we have a different period condition 
that fixes the relative height of the island with respect to other components 
of dVt. 

Let dVLi be the component of dVL around which our island forms. Let 
BVLq be a different component of dVL which respect to which the height of the 
island will be measured. Equation ([5]) implies that 

dh = 53a 



5 

TT 



and hence 



1 

height {dni) - height (a^o) = -^5 / 

Jdn^ 



dno 

a . 



where the integration is along a path connecting the two boundaries as in 
Figure [51 

For contour as in Figure [5l at the endpoints we have z = and w G M\ 0. 
In particular, the forms log z d log w and log ( d log u are integrable near the 
corresponding points of P and Q. The same computation as in proof of 
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Proposition H] shows that 

height (dfli) — height (dflo) = — Q= / logCc?logc<j + const . (19) 

cvr JdQo 

The path of the integration in fll9l) connects the points of intersection of two 
different real ovals of Q with the line z = 0. We will call it the height of one 
oval with respect to another. 

In Section I3.1.4[ we will show that the logarithmic areas of compact ovals 
of Q, the heights of the noncompact ovals, and the points of intersection of 
Q with the coordinate axes can be chosen as local coordinates on the space 
of nodal plane curves of given degree and genus. This suggests constructing 
the required curve Q by deformation from the equally weighted, no islands 
case treated in the present paper. We hope to address this in a future paper. 



3 Proof of theorem 2 



3.1 Winding curves and cloud curves 
3.1.1 Winding curves 

As we will see, the rational curve Q in the statement of Theorem [2] is not 
just any rational plane curve, it has some very special properties. We begin 
the proof of Theorem [2] with a discussion of these special properties. 
We say that a degree d real algebraic curve C C MP^ is winding if 

(1) it intersects every line L C MP^ in at least d — 2 points counting 
multiplicity, and 

(2) there exists a point Pq G MP^ \ C called the center, such that every line 
through po intersects C in d points. 

The center po in the above definition is obviously not unique: any point 
sufficiently close to Pq can serve as a center instead. The existence of a center 
implies that every winding curve is hyperbolic, that is, it lies in the closure 
of the set of curves having [d/2] nested ovals encircling the center. Much is 
known about hyperbolic curves, see e.g. j22l [23] . 

Winding curves are almost never smooth. They typically have nodes, 
but may also have ordinary triple, quadruple etc., points as well as tacnodes 
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(points where two branches are tangent). This, however, exhausts the hst of 
possible singularities in view of the following 

Proposition 5. All singularities of a winding curve C are real. Every branch 
of C through a singularity is real, smooth, and has contact of order < 3 with 
its tangent, that is, it has at most ordinary flexes. The only double tangents 
of a winding curve are ordinary tacnodes with two branches on the opposite 
sides of the common tangent. 

Proof. Let p be a complex singularity of C. The line joining p with the com- 
plex conjugate point p is real and meets C in at least 4 complex points count- 
ing multiplicity — contradicting the definition. Similarly, a complex branch 
through a singularity yields complex intersection points with any line that 
just misses the singularity, in contradiction with part (2) of the definition. 

Now suppose the origin is a real singularity of C and let a branch of C 
be parameterized by 

X = r , y = t^ + ..., 

where a < b. The line y = meets this branch with multiplicity b whereas 
the nearby line y = e meets it in 1 real point if b is odd and may miss it 
entirely if b is even. It follows that b < 3. 

The case (a, b) = (2, 3) of an ordinary cusp is ruled out by considering a 
line from the cusp to the center p. The order of contact with a nearby line 
through p will drop by 2. 

It is easy to see that a line that is tangent to two distinct points of C can 
be moved so that it loses 4 real points of contact. The same is true if two 
branches of C are on the same side of the common tangent, or tangent at a 
flex (in which case the multiplicity of the intersection with the tangent line 
is at least four greater than for that of certain nearby lines). □ 

3.1.2 Moduli of rational winding curves 

It follows from Proposition [5] that a winding curve Q is always immersed, 
that is, the map ^ 

/ : g ^ P2 (20) 

from the normalization of Q to the plane is an immersion. In this paper, we 
will focus on the winding curves of genus zero, in which case Q = P^. 

Clearly, if Q is irreducible, winding, and without tacnodes, then all small 
deformations of the map (1201) remain winding. In other words, such curves 
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form an open set of the real locus of the moduli space Mo(P^, d) of degree d 
maps 

/ : ^ P2 

considered up to a reparameterization of the domain. This moduli space is 
smooth of dimension 3d — 1, see [7]. This is because the projective plane P^ 
is convex, which, by definition, means that its tangent bundle is generated 
by global sections (concretely, the vector fields generating the S'L(2)-action 
span the tangent space to P^ at any point). The tangent space to Mo(P^, d) 
at the point Q may be identified with the global sections H^{Np2/q), where 
Np2/Q is the normal bundle (it is a bundle because Q is immersed). 

Suppose that the Newton polygon of Q is nondegenerate, that is, sup- 
pose that Q does not pass through the intersections of coordinate axes in 
P^. Consider the monomials in Q{z, w) = on the boundary of the Newton 
polygon. This gives 3d numbers modulo overall scale, matching the dimen- 
sion 3(i — 1 of the moduli space. These outside monomials determine the 3d 
intersections of Q with coordinate axes. In fact, these intersections should be 
properly considered as a point in {S'^C*Y, subject to one constraint, which 
precisely identifies it with the outside coefficients of Q{z,w). 

Proposition 6. The outside coefficients ofQ{z,w) are local coordinates on 
Mo(P^, d) near an immersed curve Q with a nondegenerate Newton polygon. 

Proof. Let ^ be a nonzero section of the normal bundle N-p2/Q. Let some 
branch of Q intersect the axis w = at the point [zq, 0) with multiplicity m. 
This intersection is preserved by ^ if and only if C, vanishes at (2:0, 0) with the 
same order m. The differential 

7] = f*{L^d\ogz Adlogw) (21) 

is a well-defined nonzero meromorphic differential on P^. One checks that, 
in fact, it is everywhere regular. Indeed, near {zq,0) we have 

w = a{z- zor + • • • , ^ = b{z- zor + ■ ■ ■ 

dw 

where a 7^ and dots stand for o{{z — -Zq)™)- Thus 

b ^ 

rj = az + . . . 

az 

and we have produced a nonzero regular differential on a genus curve, a 
contradiction. □ 
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Three types of degeneration can happen to Q in codimension one, see e.g. 
[S]. It can develop either: 

(i) a cusp, or 

(ii) a tacnode, or 

(iii) an extra node, and, hence, become reducible. 

All these degenerations mark the boundary of the winding locus in Mo(P^, d). 
By Proposition a cusp can develop only when a loop of Q ties tight around 
the center p, as, for example, for 

Q = + p=(0,0), e^O. 

In this case, the frozen boundary, which is the dual curve of Q, escapes to 
infinity and so this case will not concern us in the present paper. 

Curves with a tacnode or an extra node form codimension 1 boundary 
strata of the winding locus. A real tacnode can deform to a pair of real 
nodes or a pair of complex nodes; only in the first case does the curve re- 
main winding. Similarly, only one of the two smoothings of the new node is 
winding. 

3.1.3 Cloud curves 

The dual curve of a winding curve C separates the dual projective plane 
into the regions formed by those lines that intersect C m d and d — 2 points. 
We will call them the exterior and the interior of , respectively. The line 
corresponding to the pencil of lines though the center p lies entirely in the 
exterior of . By making it the line at infinity, the curve is placed into 
the affine plane M^. 

Since C has no cusps, has no inflection points and, hence, is locally 
convex except for cusps. The cusps point into the interior of C^, which 
makes the interior of resemble in shape a cloud, see Figure [161 This is 
why we call the dual of a winding curve a cloud curve. Node that a cloud 
curve has no real nodes other than tacnodes, for they would correspond to 
double tangents to C. 

By construction, a cloud curve has a unique, up to complex conjugation, 
complex (non-real) tangent through any point in its interior. 
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3.1.4 Higher genus 

Suppose that Q is nodal curve of genus (7 > of the kind considered in Section 
12.51 not necessarily winding. The space of nodal plane curves of given degree 
d and genus g is well known to be smooth of dimension 3(i — 1 + (7, see [8]. 
As in Section 13.1.21 above, the intersections with coordinate axes give 3(i — 1 
degrees of freedom. The remaining g dimensions come precisely from the g- 
dimensional space of holomorphic differentials as in the proof of Proposition 

m 

Proposition 7. The points of intersection with coordinate axes, the loga- 
rithmic areas of compact ovals, and the heights of the noncompact ovals, can 
he chosen as local coordinates at Q on the moduli space of genus g degree d 
nodal plane curves. 

Proof. Consider the tangent space at Q to the space of all genus g curves 
with the same outside coefficients as Q. The proof of Proposition [6] identifies 
this tangent space with the space of regular differentials rj on the curve Q. 
Moreover, by fl2T|) . the derivative of the logarithmic area (resp. height) of an 
oval Oi of Q in the direction of 77 is given by the real (resp. imaginary) part 
of / ?7 along the corresponding path /3j in Figure O 

The contours corresponding to the compact ovals of Q are closed and 
invariant under complex conjugation. Thus the corresponding periods are 
automatically real. As to the contours corresponding to the noncompact 
ovals, we can make them closed by taking their union with the complex 
conjugate contour (the dashed line in Figure E])- This makes the contour anti- 
invariant with respect to complex conjugation and thus picks out precisely 
the imaginary part of the original integral. Finally, since the contours in 
Figure [5] are disjoint, the period matrix rjj is nondegenerate. □ 

3.2 Inscribing cloud curves in polygons 

Let is a polygon in the plane, not necessarily convex, formed by 3d segments 
with slopes 0, 1, 00, cyclically repeated as we follow the boundary in the 
counterclockwise direction. Here d = 2,3, . . . . For any constant c, the image 
of Q under the map 

£xp:(a;,y)^(e^^e^^^) 

is a polygon in first quadrant M>q. The sides of £xp(n) are formed by lines 
passing through the vertices of M>q C MP^. 
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We will say that a real algebraic curve R is inscribed in the polygon 
£xp(f2) if i? is a simple closed curve contained in Q and tangent to lines 
forming the sides of fl in their natural order. Note that near a nonconvex 
corner, an inscribed curve may be tangent to the line forming a side of Q but 
not to the side itself. 

Recall that, by definition, the class of a curve is the degree of the dual 
curve. Our main result in this section is the following 

Theorem 3. Let Q be a polygon as above. A class d rational cloud curve can 
be inscribed in £xp(n) if and only if Q is feasible. If it exists, the inscribed 
curve R is unique. 

The tangent vector to R rotates d times as we go once around the curve. 
It follows that R has d — 2 real cusps, or one cusp per each non-convex corner 
of Q. Pliicker formulas, see e.g. page 280 of [4J, imply that the degree of R 
equals 2d — 2. In absence of tacnodes, in addition to the d — 2 real cusps, the 
curve R has 2d — 4 complex cusps and 2{d — 2){d — 3) complex nodes. 

The curve R being tangent to 3d lines imposes 3d incidence conditions 
on the degree d dual curve Q = R^ . As discussed earlier, one of these con- 
ditions is redundant. In fact, this redundancy translates geometrically into 
the condition that Vt is the projection of a closed contour in M?. To see this, 
note that the condition of being a closed contour is equivalent to the condi- 
tion that, as the boundary of VL is traversed, the total signed displacement 
along each of the three edge directions is zero. If the vertical edges are at 
coordinates x = xi, . . . ,Xd, the horizontal edges are at ?/ = y2, ■ ■ ■ ,yd and 
the slope-1 edges are ed, x — y = Zi, Z2, ■ ■ ■ , Zd, then the displacements along 
horizontal edges are yi + Zi — Xi,y2 + Z2 — X2, ■ ■ ■ ,yd + Zd — Xd- To be closed the 
sum of these displacements must be zero. The fact that the displacements 
in the other two directions are zero gives an identical relation. The equation 



is the logarithm of the condition that the intersection points of a plane curve 
with the three coordinate axes has product 1. 

Also note that while imposing 3d — 1 incidence conditions on a ratio- 
nal plane curve of degree d gives generically a finite set of possibilities, the 
number of possibilities grows very rapidly. For example, there are 26312976 
rational curves of degree 6 through 17 general points in the plane (all of these 



d 




1=1 
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numbers were first determined by Kontsevich, see e.g. [7]). This shows that 
solving the incidence equation requires some care. In fact, our existence proof 
can be turned into a practical numerical homotopy procedure for finding the 
unique inscribed curve. 

3.3 Proof of Theorem [3] 
3.3.1 Strategy 

Our strategy will be deformation to c — > ±cx) limit, in which tropical alge- 
braic geometry takes over. This is essentially the patchworking method of 
O. Viro, see e.g. [9J. A closely related approach was employed, for example, 
in [16] to construct curves of given degree and genus passing through given 
points in the plane. 

From the point of view of random surfaces, letting c — > ±oo means im- 
posing an extreme volume constraint. As c ^ +oo, the limit shape becomes 
the unique piecewise linear surface that minimizes the volume enclosed, see 
an example in Figure O 




Figure 6: Surface enclosing the minimal volume 

In the c —>■ ±oo limit, the inscribed curve collapses to (a part of) the 
corner locus (locus where the surface is not smooth) of this piecewise linear 
surface, an example of which is plotted in yellow in Figure El The relevant 
parts of the corner locus are convex edges for volume maximizers and concave 
edges for minimizers (as seen from above the graph). 

Our goal is to reverse this degeneration, that is, starting with the concave 
corner locus of the volume minimizing surface, we construct a certain tropical 
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curve Qt, the would-be dual of the inscribed curve. Then we show that Qt 
can be deformed to a rational degree d winding curve Q{x,y; c) = whose 
dual curve R{x,y;c) = is inscribed in £xp(n) for all sufficiently large c. 
Finally, we prove that these curves can be further deformed to any finite 
value of c. 

In order for the very first step to work, we need to assume that Q is 
feasible, otherwise there is no volume maximizing surface to begin with. We 
further assume that Q is generic, that is, there are no accidental relations 
between its side lengths. In particular, for generic fl, the triple points of 
the concave corner locus of the volume minimizers are local minima, see the 
illustration in Figure [71 Moreover, genericity implies that the maximizing 
surface lies strictly above the minimizing surface, except at the boundary. 
This implies that there are no "taut" paths, that is, paths in the interior 
which are present in every surface. See Figure [8] for an example of a taut 
path. 





Figure 7: A nongeneric polygon fl and a nearby generic one. 



3.3.2 The tropical curve 

The construction of the tropical curve is very simple: take the concave corner 
locus of the volume minimizer and extend univalent and 2-valent vertices to 
3-valent ones as in Figure [9|, where the newly added rays are shown in blue 
(thinner lines). 

The tropical curve Qt thus obtained is a genus zero degree d tropical 
hyperbolic curve, see [22] where properties of such curves are discussed in 
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Figure 8: A polygon with a taut path (the vertical bisector). In a tiling, no 
tile crosses the taut path. 



Figure 9: The tropical curve associated to the surface in Figure [61 
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detail. 

Consider the connected components of the complement of Qt, which we 
will call chambers. There are {^^^ chambers and they naturally correspond 
to lattice points in the triangle with vertices (0,0), ((i, 0) and (0,(i). The 
meaning of this triangle for us is that it is the Newton polygon of the degree 
d curve y; c). 

The dual of the graph Qt is, combinatorially, a tessellation of this triangle 
with unit size triangles and rhombi (recall that we assume Q, is generic). Each 
rhombus in the tessellation will correspond to a node of Q^x., y; c) = 0; there 
will be (^-^) of such. 

Note that triple points of Qt come in two different orientations, and triple 
points of one of the two orientations correspond bijectively to local maxima 
of the concave corner locus. 



3.3.3 Construction of Q 

We will first construct an approximation 

Q°(x,y;c)= Yl ^S^V' (22) 

i+j<d 

to the true winding curve Q. 

Recall that the monomials of Q correspond to chambers of Qt- For exam- 
ple, monomials that differ by a factor of x correspond to adjacent chambers 
separated by a vertical line (which sometimes can have zero length). Let 
X — aij be the line separating the chambers corresponding to the monomials 
x^y^ and x^'^^y^ . Then we require that 



Similarly, for monomials that differ by factor of we set 



if the line separating the chambers is y — fiij. For monomials which share a 
diagonal edge x — y = jij we have 
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Since at a trivalent vertex we have 



f^ij Oiij 

this rule defines the curve (5° = uniquely and unambiguously up to a 
common factor. The irrelevant common factor may be chosen so that all 
coefficients of are real and positive. 

By construction, as c oo the amoeba of the curve Q^, scaled by c, 
converges to the tropical curve Qt- In particular, intersects the coordinate 
axes approximately at the points 

(l,-e-^^%0), (23) 

where 

x = au y = Pi, x-y = -fi, i = l,...,d. 

are the lines forming the 3d sides of Q By a small c-dependent change of the 
original curve Qt, we can make the curve Q'^ defined by ( !22|) pass precisely 
through the points (!23|) . 

A small deformation of a quadruple point of Qt (into a pair of triple 
points connected in one of two ways) leads to the change of the topology of 
the real locus of Q^ for all large c. There is a c-dependent deformation of Qt 
for which every quadruple point becomes a real node, see [16]. Thus, for all 
c ^ 0, we have constructed a genus zero curve passing through the points 
f l25]) . This is our curve Q{x,y; c). 



3.3.4 Intersections with lines 

Our next step is to prove that the curve Q{x, y; c) is winding. This involves 
understanding its real locus. We know that the curves 

Q(±e^^±e"*';c) =0 (24) 

approach Qt as c ^ oo. One of these four curves, the one that corresponds 
to the (+, +) choice of signs is empty because Q has only positive coefficients. 
The shape of the other three curves can be inferred from simplest example 
of a conic, shown in Figure [TOl The colors in Figures [TOl [12] and [13] are 
explained in the following table: 



quadrant 



color 



[+,+) (+,-) (-,+) 



dotted black dashed gray 
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The essential feature seen in Figure [TU] is that the black and dashed curves 
cross over near the middle of the corresponding edge of the tropical curve. 
(Note that for clarity in some of the figures we use rectilinear coordinates and 
in others, coordinates with the 3-fold symmetry. Hopefully these coordinate 
changes will not confuse the reader). 



Figure 10: The curves Q(±e^, ±e^) = for Q{z, w) = l + z + w + e ^zw. 

This same feature holds in general since for c sufficiently large, the struc- 
ture of the curves near a compact edge of Qt only depends on the coefficients 
of Q corresponding to the four chambers around that edge: the other coeffi- 
cients are exponentially small compared to these. It follows that the curves 
f l24p . in our running example, look like the curves in Figure [TTl 

Now we turn to intersecting the curves Q{x, y; c) with lines. There are two 
kinds of lines, those that intersect the positive quadrant and those that don't. 
The latter ones can be defined by an equation with positive coefficients. The 
difference between the two kinds of lines is illustrated in Figure [121 As we 
will see only lines that miss the positive quadrant have a chance to intersect 
Q{x, y; c) in fewer than d points. 

A basic property of any generic tropical curve, and of Qt in particular, 
is that it intersects a generic tropical line in d points, where d is the degree. 
Generic here means that the triple point of the line is not a point of the 
tropical curve. It is easy to see that in this generic situation there are still d 
points of intersection of the nearby nontropical curve Q with a nearby line: 
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Figure 11: This is how the curves (5(±e^^, ie'^^; c) = look hke for c 3> 
for the tropical curve in Figure M 
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transverse intersections are preserved under deformations. 

If a line is not generic, some of these points of intersection may be lost. 
Note that real points have to be lost in pairs, so, for example, if the triple 
point of the line lands on a (noncompact) ray of Qt, there are still d real 
points of intersection for nearby curves (this follows because the point at 
infinity has two colors for both Q and the line and so this point is still a point 
of intersection). However, if the triple point of the line lands on a compact 
segment of Qt, two points of intersection may disappear. If the triple point 
of the line lands on a triple point of Qt of the opposite orientation, two 
points of intersection (at most) may disappear. Finally, if the triple point of 
the line is a triple point of Qt, with the same orientation, and all incident 
segments of Qt are compact, then as many as 4 real points of intersection 
may be lost (Figure [T3l). 




Figure 13: How to lose four intersections. The Q curve here is in thinner 
colors, the line in thicker colors. 

A simple but crucial observation is that, by construction, the curve Qt 
has no triple points of the same orientation as the tropical line and incident 
to 3 compact edges. This is because, as discussed earlier and illustrated in 
Figure [7l triple points of compact edges correspond to local minima of the 
volume minimizer for generic Q. 

Furthermore, it is elementary to check that if the line intersects the pos- 
itive quadrant then it has d points of intersection with Q. Similarly, any 
tropical line with a triple point sufficiently far away from the amoeba of Qt 
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meets Q in d points. Thus the curve Q is winding. 

The dual cloud curve R = Q^, which, recall, is convex and non-sef- 
intersecting, is approximately the yellow convex corner locus from Figure [61 
This is because the tropical lines whose triple points lie on a compact edge 
of Qt are exactly those which lead to fewer than d points of intersection. 
The curve R is inscribed in Q because, if it were protruding it would have 
more, say, vertical tangents than we constructed. This is impossible since 
more than d vertical tangents means more than d points of intersection of Q 
with a certain line. 

3.3.5 Deforming the winding curve 

We will now show that the inscribed curve that we constructed for c ^ 
can be deformed to any given value of c. Changing c produces a particular 
1-parameter family of polygons into which a cloud curve is to be inscribed. 
It is more convenient, in fact, to consider the general variation of the polygon 
£xp((]). 

In the set of all polygons Q with 3d disjoint edges (that is, no part of dfl 
is traversed twice) consider the set of all feasible ones. Denote the interior of 
this feasible set by T. Each connected component of T is naturally an open 
set of R'^'^"^. Let Tq C T be the set of those Q for which there exists an 
irreducible cloud curve R without tacnodes inscribed in £xp(i7) with c = 1. 
Our goal is to prove that Tq = T. 

We know that the intersection of Tq with each connected component of 
T is nonempty. In fact, if Q is generic, then cQ, c ^ 0, is contained in Tq. 

Observe that Tq is open. This follows from Proposition because chang- 
ing Q amounts to moving the points of intersection of Q with the coordinate 
axes. 

We claim that Tq is also closed. Indeed, let be a limit point of Tq. Since 
the set of winding curves is closed and R remains bounded, there will be an 
inscribed cloud curve R in £xp(n), but it might be reducible or tacnodal, 
or both. Since the map from R to the circumscribed polygon is continuous, 
it is enough to show that the codimension 1 degenerations cannot happen, 
that is, it is impossible to acquire exactly one tacnode or one extra node. 

If the curve R develops a tacnode this means that a part of the frozen 
region is swallowed by the liquid region as in Figure O This is impossible 
because, for example, there must be a point of tangency with the boundary 
between any two cusps of R. 
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Figure 14: A cloud curve develops a tacnode 



The other potential degeneration is when Q develops an extra node (and, 
hence, becomes reducible). First suppose that neither component of the 
degenerate curve is a line. In this case, the dual curve R becomes the union 
of two cloud curves and a double line tangent to both of them. Nearby 
cloud curves look like a union of the two disjoint cloud curves connected by 
a thin elliptical piece along a double tangent, see illustration in Figure [15] 
on the left. This, for example, is what happens if we squeeze the polygon 
in Figure [6] in the middle until the middle portion goes to zero width. This 
example, however, is not generic because the linear pieces of the limit shape 
on two sides of the collapsing ellipse have different slopes and meet along a 
line in a coordinate direction. Generically, the double tangent is not a line 
in a coordinate direction and these linear pieces (which then have the same 
slopes) merge together. This leads to a taut path in the polygon and, hence, 
cannot happen in our situation. 




Figure 15: A cloud curve becomes reducible 



The other possibility is when Q degenerates to a union of a winding curve 
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Qi and a line Q2- The line Q2 corresponds to a point p2 of the dual plane and 
the collapsing ellipse follows the tangent from p2 to Q^, this is illustrated in 
Figure [15] on the right. This degeneration can only happen if an edge of fl 
has length tending to zero. 

Note that in codimension two one can have a degeneration in which two 
cusps of the frozen boundary come together and merge in a tacnode, remi- 
niscent of the Henry Moore's Oval with Points. This corresponds in Figure 
[T5| left part, to the thin ellipse having zero size. 

Finally, there is practical side to the above deformation argument. Namely, 
it allows us to find the inscribed curve by numeric homotopy, that is, by solv- 
ing the equations using Newton's algorithm with the solution of the nearby 
problem as the starting point. An example of a practical implementation of 
this can be see in Figure [T6l 




Figure 16: Frozen boundaries (at c < 0, c = and c > 0) for the polygonal 
region in Figure [61 

3.3.6 Uniqueness 

It remains to prove that the inscribed curve is unique and that the polygon 
is feasible if the inscribed curve exists. Given the inscribed curve R, let 
denote the height function obtained from it by the procedure of Section 11.81 
This is a well-defined function with gradient in the triangle ([1]). Therefore Q 
is feasible. The inscribed curve is unique because by the results of Sections 
13.3.51 it has a unique deformation to the unique tropical inscribed curve. 
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4 Global minimality 



As before, let h* denote the height function constructed from the inscribed 
cloud curve R. 

Theorem 4. The height function h* is the unique minimum for the varia- 
tional problem (j4]). 

Proof. By strict convexity of the surface tension, there is a unique minimizer. 
Therefore we need to prove that 

m > (25) 

for any function g satisfying the boundary conditions and with gradient in 
([1]). Since the integrand in (jlj) is bounded, we can freely modify it on sets 
of arbitrarily small measure. In particular, it is enough to prove (1251) for 
functions g such that g = h* in a. neighborhood of the vertices of Vt. 

For any such g, we will construct a sequence {/i„} converging to h such 
that 

(i) £(/i„) £(^*), as n — i> oo , 

(ii) £ has a directional derivative at hn in the direction oi g — km 

(iii) this directional derivative is nonnegative for all n ^ . 

By convexity of £, the last condition implies that E{g) > 8.{hn) for all n ^ 0, 
from which (125|) follows. 

Now we proceed with the construction of the sequence {hn}- In the 
liquid region we will take hn = h*, except in a very thin strip along the 
frozen boundary. Inside the facet, hn will be a smooth function very slightly 
deviating from h* and equal to it in a fixed neighborhood of the vertices of 
Q. Note that the facets of h* are of two types: the "bottom" facets, for 
which g > h* for any allowed function g, and the "top" facets, where the 
inequality is reversed. For example, the bottom facet in Figure [1] is also a 
bottom facet according to the above definition. For concreteness, we will 
give a construction of hn on a bottom facet; the case of a top facet is the 
same, with obvious modifications. 

The construction of Section 13.31 give a family of solutions of the Euler- 
Lagrange equation ([3]) parameterized by a real constant c. Since we will 
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need all of them for the argument that follows, we will use Cq to denote the 
particular value of the Lagrange multiplier in (jl]) for which we wish to prove 
that h* is a minimizer. Note that the height at any point is a decreasing 
function of c. Indeed, the c-derivative = ^ of the height function solves 
the linearization of ([3]), which has the form 

divG- V0 = 1, 

with a positive definite matrix G = V^cr(V/?.*) in the liquid region and zero 
boundary conditions on the frozen boundary. The condition G > guaran- 
tees (j) has no local maxima and hence (f) < 0. 

Consider a horizontal bottom facet of h*. Its boundary is formed by a 
part of dQ and a part of the frozen boundary. Let us denote that part of 
the frozen boundary Fc^ and call it the frozen front. For c below cq, the 
monotonicity ^ < implies the frozen front Fc moves inside the facet. It 
sweeps the entire facet as c — * — oo. 

By our assumption on g, we can find C <^ such that g = h* outside 
Fq. The functions will smoothly interpolate between h* outside Fc and 
h* just inside in such a way that 

< \\VhJ <ess 0(n~^/2) ^ n^oo. 

between Fc and F^g. Here <ess means that the inequality is satisfied except 
on a set whose measure goes to zero as n oo. 

The derivative of £ at /i„ in the direction g — hn equals 

j (c-divVa(V/i„))((7-/i„), 

which is nonnegative for all sufficiently large n provided g > h* somewhere 
on the facet and 

div Vo-(V/i„) < c 

between Fc and Fcg. 

The functions /i„ will be obtained by patching together solutions of 

divVa(V/,) =Q, (26) 

where 

C = Co > Ci > C2 > C3 > ■ ■ ■ > = C , 
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are chosen so that all the distances between the fronts F^- go to zero as 
0{n~^) as n ^ oo. The implicit constants in 0(. ..) here and below are 
bounded in terms of C. 

For each Cj, consider the algebraic solution of fl26|) with frozen front Fc.. 
Since near the frozen boundary its gradient has a square root growth, we have 
||V/i|| = 0(^-1/2) between F^^ and F^^_^. This means that the discontinuous 
surface formed by the union of /j's between Fc^ and Fc-_-^ is nearly flat for 
n 0. To make it continuous, we shift each successive piece slightly in the 
X3 direction (recall that our facet is assumed to be horizontal) and insert 
a smooth nearly vertical strip of surface between two /j's. The width of 
each of these strips is O^n^^^"^) and so their contribution to the functional is 
negligible. There will be also a small region near the boundary of Q where 
fi is not defined. Since each Fc- is tangent to dfl, the area of that region 
is 0{n~^) and so we can extend fi there as an arbitrary smooth Lipschitz 
function. Thus we have constructed the required sequence {/i„}. □ 

5 Explicit examples 

In this section we work out some explicit solutions to ffTOl) for various cases. 

It can be checked along the lines of Section H] that the solutions we con- 
struct below are minimizers for their corresponding variational problem. 

5.1 Disconnected boundary example 

The simplest multiply-connected regions are obtained from simply-connected 
regions by removing points. In this case the inner boundary components are 
reduced to points and specifying a surface with these boundary conditions 
is equivalent to specifying the height of the surface at given points in the 
interior of the domain. 

Figure [T71 shows an example: here the outer contour is a hexagon (6 edges 
of a cube) and the inner contour is a single point, the central point in the 
figure, located at a height which is 2/3 of the way between the upper and 
lower vertices of the corresponding cube. 
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The curve Q is a cubic curve which by symmetry has the form 

+0W +azw +bz w 

+ 1 +bz +bz^ +z^. 

The coefficient b is determined by the intersection of Q with coordinate axes, 
that is, it determined by the size of the cube and the value of Lagrange 
multipher c. The coefficient a determines the height of the central point. 



Figure 17: Boxed plane partition conditioned to pass through a given point. 
The gray edges lie below the surface. 



5.2 Example with a bubble 

In this section we consider dimers on the square-octagon graph. Here P{z, w) = 
S + z+l/z+w + l/w, see [I3j. The Newton polygon is defined by |x| + < 1. 
Consider the contour C through the points 

(1,0,0), (0,1,1), (-1,0,0), (0,-1,1), 

see Figure fT8l 

We expect 4 frozen phases and hence 4 points where the frozen phases 
change along the frozen boundary. Such points correspond to points of P at 
the toric infinity, that is, tips of tentacles of the amoeba of P. We, therefore, 
expect Q to have the same Newton polygon as P. From symmetry, we 
conclude that 

Q{z,w) = b + z + 1 / z + a{w + 1 / w) . 
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The point {z,w) = (00,00) on P occurs somewhere along the upper right 
boundary edge (the hne x + y = 1) since this is the point of separation of the 
two adjacent frozen phases. Ehminating w using ffTOl) . substituting y = 1 — x, 
and expanding near z = 00 gives 

where (xq, 1 — xq) is the point of intersection of the frozen boundary with the 
hne X + y = 1. This leaves us with one variable, xq, which is determined by 
the area constraint from Proposition HI 

While it is a non-trivial computation to compute Xq in general, in the 
case c = we have an extra symmetry and expect x = 1/2. This leads to a 
degree 8 equation 

729 - 13608a;2 - 22896a;^ + 64000x^ + 102400a;^ - 13608?/^+ 
412992a;\2 - 1104000x^?/2 + 870400x^^2 - 22896?/^ - 1104000x\^+ 

2054400xV + 64000?/^ + 870400x^7/^ + 102400?/^ = (27) 

for the frozen boundary in the unconstrained volume case-see Figure [TSl 
This was called the "octic circle" by Cohn and Pemantle [2]. 

This case can be obtained more simply using IQ: there let Qo{z,w) = 
-3/2. We solve 

P{z,w) = 
3 

-- = xzPz + ywPu, 

for z,w as a. function of x, y. The double root of these equations occurs along 
fl27|) . as one immediately checks using resultants. 

5.3 Crystal corner boundary conditions 

In this section we consider minimizers defined on the whole plane and satisfy- 
ing piece-wise linear boundary conditions at infinity. They can be interpreted 
as dissolution shapes of an infinitely large crystal near one of its corners. 

5.3.1 Minkowski sums of Ronkin functions 

Recall that the Minkowski sum of two sets A,B E M" is 
A + B = {x + y\xeA, y e B}. 
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Figure 18: frozen boundaries for the square-octagon "fortress" example. 



If A and B arc convex, so is their sum. 

For a curve Q{z,w) let Aq be the set of points in lying on or above 
the graph of the Ronkin function JIq of Q. This is a convex set. 

A general class of solutions to the volume-constrained minimization prob- 
lem for an arbitrary spectral curve P can be constructed as follows. 

Theorem 5. Let P' he a Harnack curve with the same Newton polygon as 
P such that the areas of the corresponding (compact) facets of !Jlp and !Jlp/ 
are equal. Then the boundary of the Minkowski sum of the sets Ap and Api 
is a minimizer. 

Proof. It is a simple geometric property of the Minkowski sum that a facet on 
the boundary of Ap + Api is the Minkowski sum of a facet of and a facet 
of Api having the same slopes, where one of these facets may be reduced to 
a point. Any point p on the boundary oi Ap + Api but not on a facet is the 
sum of unique points on the boundary oi Ap and Api with tangent planes in 
the same direction. 

Given a point (xi, yi, ^Kp(xi, yi)) on the graph of 3?p, which is not in 
the interior of a facet of Dip, there is a corresponding point (2:1, Wi) on the 
curve P, with log 1 2:1 1 = Xi and loglwil = 7/1, and such that the arguments 
of z\ and wi are linearly related to the slope of "Jip at {x\,y\). The point 
(2:2, 1/2, 3?p/(a;2, ^2)) with the same slope on the graph of IRp/ (and if this 
point is on the boundary of a facet, it must also have the same slope along 
the facet boundary as the corresponding point of Jlp) yields a corresponding 
point (2:2,^2) lying on the curve P' and with the same arguments. 
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Let Q{z,w) = P'{l/z,l/w). Construct the minimizer using the curve 
Q in (fTU]) . Since Qi^e''^^ Zi, e~'^ywi) = 0, we have (-22,^^2) = {e'^^ / zi,e'^y /wi). 
Therefore 

(2:2,1/2) = (I0g|z2|,l0g|w2|) = 

(cx-log|2;i|,cy-log|wi|) = c - {x,y)- 

and so {x,y) = + (x2,l/2)). 

The scaled Minkowski sum is the unique shape whose tangent plane above 
{cx,cy) has the same slope as the tangent planes at {xi,yi) and {x2,y2)- 
Therefore it is the solution to the minimization problem flTUl) . except that it 
may not satisfy the area constraints of Proposition HI Those say precisely 
that the facet areas are equal. □ 

In \TT\ we proved that, for fixed boundary behavior, there is a unique 
Harnack curve with facets of given areas. So we can always find P' with the 
appropriate facet areas. 

Note that if the facet areas of P and P' do not match, we still get a C^/^ 
function that solves the Euler-Lagrange equation in the liquid region and is 
locally a minimizer inside each facet. It is not a global minimizer since raising 
or lowering a facet as a whole will reduce the total surface tension. From the 
point of view of Glauber dynamics on random surfaces, moving the facet as 
a whole is an extremely slow process, so this solution of the Euler-Lagrange 
equation may be described as metastable. 

Note also the symmetry between P and P': Ap + Api is a minimizer for 
the variational problems for both P and P'. 

There exist intriguing connections, first noticed in [121, between the ran- 
dom surfaces studied in this paper and topological sting theory, see e.g. [18] 
for a review aimed at mathematicians. In this context, the Minkowski sum 
construction has a natural interpretation which will be explained in [T^ . 

5.3.2 Wulff construction 

Taking P' = P in Theorem O we have Ap = Api and the facets trivially have 
the same area. So limit shape is a copy of the Ronkin function of P. This 
shows that the Ronkin function itself is a solution to the volume-constrained 
surface-tension minimization problem. Indeed, this is the classical Wulff 
construction, since the Ronkin function is the Legendre dual of the surface 
tension [13] . 
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5.3.3 Higher-degree solutions 

To combine P with higher-degree curves in Theorem [5|, one can take covers 
of P as described in [13]: the curves 

have the same Ronkin function up to scale as P, and N{Pn) = nN{P), where 
N[P) stands for the Newton polygon of P. 

Let Q be a Harnack curve with N{Q) = N{Pn), and having facets of the 
same area as Pn- Then Ap^ + Aq is a scaled minimizer for P„, and since 
ap^{ns,nt) = ap{s,t), Ap^ + Aq is also a scaled minimizer for P. 

This allows us to combine P with arbitrarily high-degree Harnack curves, 
and in fact by taking limits one can construct explicit non-algebraic mini- 
mizers. 



5.3.4 Rational curves 

A particular case of Theorem O is when P has genus zero — then Q must also 
have genus zero. Suppose P{z, w) = z + w + 1. Then P„ can be parameterized 
via z I— > (2" : (—1 — z)" : 1). A degree-n rational Harnack curve Q has a 
parametrization 



n 

\ A Y[{u -ai):B - h) : C - c,, 

i=l i=l i=l 

with A,B,C G M and where the real numbers oi, . . . , a„, 61, ... , bn, ci, . . . , c„ 
occur in cyclic order around M U 00 (that is, all the a's come first, then the 
6's and then the c's), see [TT]. 



Then (ITUI) gives 



„ n 

A^^ u- a, 

i=l ' 



„ n 



B u-bi 

The constants c. A, B, C can be incorporated into an overall scale and 
translation. Suppose that as n ^ 00, the measures 
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have weak limits that we will denote by (i/ij, i = 1,2, 3. Taking logarithms 
and limit, we find that the triangle discussed in Section takes the form 

J2 e""""'^"' exp J \og{u - t)dni{t)^ = , (28) 
for certain constants Cq, . . . , C3. 
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